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SUMMARY 


The indiclal lift and pitching-moment coefficients are derived for 
flat— plate triangular wings traveling at supersonic speeds. The coeffi- 
cients are determined for angle— of-attack distributions corresponding to 
sinking wings and to pitching wings. The wing with supersonic edges is 
completely analyzed, and the wing with subsonic edges is partially 
analyzed, the solution in this case being completed for very narrow 
wings by an application of slender wing theory. In the case of the 
supersonic edges, a comparison is made with known two-dimensional 
results and also with the results for the same triangular wing in 
reversed flow. 


INTRODUCTION 


The wing of triangular plan form has received considerable atten- 
tion in the steady— state theory of three-dimensional wings in a super- 
sonic stream. The purpose of the present report is to determine the 
aerodynamic characteristics of a triangular wing in supersonic unsteady 
motion. 

There are several simple types of unsteady motion on which the 
analysis can be based. The so-called indicial motion, in which the 
velocity undergoes a discontinuous change at t*=0, will be considered 
here. (See also references 1 and 2 .) It is possible to conceive the 
physical situation in two slightly different ways. For one, it can be 
supposed that the wing has been traveling at the constant velocity 7 q 
for an infinitely long time and then, at t*=0, starts suddenly to sink 
without pitching motion (or to pitch without sinking) while maintaining 
the forward velocity Vq, On the other hand, the wing may be considered 
to be at rest in still air until at t'=0 it starts suddenly either to 
sink or to pitch and, at the same instant, attains the forward velocity 
Vq. The latter physical picture will be used in this report. Problems 
of unsteady motion can also be approached with the initial assumption 
that the velocity potential depends harmonically on the time, (See 
reference 3.) These two approaches are quite compatible in that they 
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can te related through the use of superposition methods (iXihanBl's 
integral, Fourier's integral) of the operational calculus. 

The question of whether the coordinate system should move with the 
wing or remain fixed is also of some importance. (See reference h for 
a discussion. ) The latter alternative, that is, where the wing moves 
away from the coordinate system, was chosen for this report because the 
velocity potential in this case satisfies the wave equation 

'Pxx + 'Pyy 'Pff = ® 

“o 

where x,y, z are Cartesian coordinates, t* is time, and a© is the 
speed of sound in the undisturbed medium. The fact that the equation 
has this form is helpful in establishing analogs between steady and 
nonsteady motions, and these analogs are of considerable help in the 
solution of certain problems. 

The boundary conditions to be considered correspond to the problem 
of the flat sinking wing (angle— of— attack distribution uniform over the 
plan form) and to the flat pitching wing (angle of attack varies 
linearly with chordwise distance) in indicial motion. First, the load- 
ing on a flat triangular wing with supersonic edges undergoing an 
indicial sinking motion is determined. Then a simplified method is 
developed whereby total lift and pitching-moment coefficient for the 
wing with supersonic edges may be obtained. These quantities are 
determined as functions of time, for both sinking and pitching wings. 

Lastly, the triangular wing with subsonic edges is partially 
analyzed, and an approximate method for very slender wings is used to 
complete the determination of loading. The analog method, mentioned 
previously, is here of great value. Lift and pitching-monent coeffi- 
cients for the sinking and pitching slender triangular wing are 
determined. 


LIST OF IMPOETANT SYMBCLS 


ao speed of sound in the free stream 
Cq root chord of triangular wing 

Cl lift coefficient ( . 

\ho Vo^S 




■) 


’a 


indicial lift coefficient due to angle-of-attack change (without 
pitching) 
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CL(j' indicial lift coefficient due to pitching on a wing rotating about 


its leading edge or apex 


■ 


8Cl 




V- ' . 

*0 fl_ 


0 = 0 / 


Cm pitching-moment coefficient, positive when trailing edge tends to 
sink relative to leading edge 

indicial pitching-momenb coefficient due to angle-of-attack change 
(without pitching) measured about the leading edge or apex 

8Cti 


^ ^ a=0^ 


indicial pitching-moment coefficient due to pitching measured about 

\ 


the leading edge or apex ( 


8 Cm* 

. Cq 0 

^ Vo 5=0. 


m 


cotangent of sweep angle (cot A) 

Mq free-stream Mach number ( ) 

\ Sq y 


Mq 

2as 




<iy 


2 0 s^ 


<lo 


'lo 


loading coefficient (pressure on the lower surface minus pressure 
on the upper surface divided by free— stream dynamic pressure) 


free-stream dynamic pressure 
q dimensionless rate of pitching 


(^) 


S 

S 

t* 

t 


local semi span of wing 

wing area 

time 


«o^’ 
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■to 

U,T,W 

Vo 

x,y,z 

a 

P 

Pe 

e 

0 

A 

T 

<P 

A(f> 


n 

u 


t 

Co 

perturbation velocity components in the x^y^z directions, 
respectively ' 

free-etream velocity 

Cartesian coordinates 

distance of center of pressure from wing apex 

angle of attack (angle between fli^t path and plane of wing), 
radians 

yii-Mo^i 

1 

m Mq 4 */l+vF 

wing qngift of pitch, relative to initial attitude, positive when 
trailing edge lies below leading edge 

wing rate of pitch, positive when trailing edge is sinking 

relative to leading edge 

angle of sweep of leading edge, positive for sweepback 

free— stream density 

t 
s 

/ Vot’ 

chord lengths traveled I — • or Moto 

\ °o 

perturbation velocity potential 

Jump in potential across the z=0 plane 
[<P(x,y,0‘'’) - cp(x,y,0 )] 

Subscripts 

component taken normal to the leading edge 

positive side of the z=0 plane, or upper surface of a wing 
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THE 


PROBLEM 
OF ITS 


AND THE NATURE 
SOLUTION 


Consider a wing situated in still air and a Cartesian coordinate 
system with origin on the leading edge at the point of wing symmetry. 

At time equal to zero, the wing staarts impulsively to move in a straight 
line with constant velocity away from the coordinate system which 
remains fixed relative to the still air at infinity. The load distri- 
bution on a wing undergoing such a motion is called the indicial loading. 
Similarly, the forces and moments which are based on this loading are 
given the adjective indicial. The partial differential equation that 
is satisfied by the velocity potential 9 for such a motion can be 
written 

nt -'Pn-'l’yy-'Pzz - ® (1) 

which is the normalized form of the wave equation. In equation (l), 
x,y, and z are distances; x measured chordwlse, y spanwise, and z 
vertically, and t is equal to aQt' where aQ is the free— stream 
si>eed of sound and t* is time. ^ 

The boundary conditions to which equation (l) is subject are 
dependent on the wing shape and motion. Adopting the assumptions of 
thin-airfoil theory, which eu*e consistent with the assiui 5 )tions already 
used in obtaining equation (l), it can be assumed that the slope of the 
wing surface, in the direction of motion, at any place and time is 
given by the ratio of tl» vertical velocity component in the z=0 plane 
to the wing’s forward velocity component. For a flat plate, then, the 
following conditions are to be satisfied;- 

1. The vertical velocity 9^ is a llneeir function of x, 
the coefficients of which depend upon the angle of attack 
rate of pitch, over the portion of the xy plane occupied by 
the wing at any given time. 

2. No perturbations exist at infinity. 

3 . There are no discontinuities in the velocity potential 
except over the region occupied by the wing and its vortex wake. 

The problem has now been expressed as one of finding, for pre- 
scribed boundary conditions, a solution to the .wave equation. It is 
often desirable to express the solution in terms of the loading coeffi- 
cient rather than the potential function or velocity components. This 
coefficient can be written in its linearized form as 
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^ ^ 2 8ACp , 

<lo VqMq "SF 

The solution of the problem without further restrictions has, in 
general, not been obtained, although it can be shown that the boundary 
conditions for any wing in unsteady motion can be satisfied by a 
suitable superposition of sources and doublets (reference 4). However, 
the solution for wings with all supersonic edges can be written as a 
double integral of sources having intensities determined by the local 
slope of the wing. Hence, for a flat surface at constant angle of 
attack and not pitching. 


9(x,y) 



i dxidyi 


(3) 


where is the outline of the region of sources which can, at a 

given time, affect the point x,y at which the potential is being 

determined, r equals V^(x— Xi)^ + (y-yi)^, and a is the angle of 
attack of the wing. The area Sg^ has been termed the acoustic plan 
form eind a discussion of its significance is given in reference 4. 


The. solution for the triangular wing with subsonic edges can be 
obtained in certain regions, but in others the problem reverts to the 
solution of a double integral equation involving time and the two sur- 
face dimensions of the wing. If the triangular wing is slender, an 
approximate method for finding the pressure over the entire wing can be 
used. This method is to neglect the streamwise velocity gradients in 
comparison with the gradients in the plane normal to the free stream 
and also in comparison with the term 9.(;t* There results for the 
partial differential equation governing the flow field (equation (l)) 
the wave equation of one lower dimension, namely, 

■ftt - - V = 0 C*) 

As will be developed later, the boundary conditions become the same as 
for a rectangular flat plate of very low aspect ratio inclined at angle 
of attack to a free stream with a Mach number equal tO\/^ This 
analogy with the steady-state lifting— surface problem is useful since 
solutions to the latter problem have been obtained. 

The indie ial lift and pitching-moment coefficients will be deter- 
mined for wings with two different vertical velocity distributions. 

The first is the case in which the boundary condition is 
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vhere a is the angle of attack, and is constant over the plan 

form. This case is referred to as that of the sinking wing. The other 
boundary-value problem considered is the one corresponding to a flat 

pitching at a constant rate about its leading edge or apex. The 
boundary conditions have the form 

(x+Mot) 

where 0 is the constant rate of pitch and (x + Is the distance 

aft the leading edge or apex of the wing. The pitching velocity 0 
is considered positive when the trailing edge of the wing sinks relative 
to the axis of rotation. 

The difference between a 
and 0 is illustrated in the 
sketch, where for clearness the 
motion is shown to be oscillatory 
rather than indie lal. The angle 
of attack a is the angle 
between the flat wing surface and 
the line tangent to the fl 1 gh t 
path of the leading edge or apex 
of the wing. The angle 0 is 
the angle between the flat wing 
surface and the horizontal (see 
part (a) of the sketch). Part (b) 
of the sketch shows a wing under- 
going a sinusoidal angle-of-attack 
variation with a zero angle of 
pitch throu^out the motion. 

Part (c) shows a wing undergoing 
a sinusoidal angle— of— pitch varia- 
tion taken about the leading edge, 
with the angle of attack remain- 
ing zero. 

Although the lift and 
pitching-^noment coefficients are 
given only for the two types of 
motion separately, it is possible, 
according to the linear theory, 
to combine them to simulate an 

arbitrary indlcial maneuver consisting of both sinking and pitching 
motions. For such a case, the boundary condition becomes 

Wu = — -0 (x + Mot) 


Direction of wing motion 
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If, furthermore, the desired maneuver is not indicial, the lift and 
pitching-moment coefficients may he determined by use of Duhanel’s 
integral : 

(t'-T') + C^*(f-T*k( T*)] dT« 

j n t ’ 

Cm = )a( T') + Cj^’(f-T') q( T*)] dT* 

where cx(t*) and q.(t*) are the arbitrary motions (• ■ ty %. / Sx'" 

CLqS and Cmg^’ are the indicial aerodynamic coefficients, and the primes 

on the coefficients indicate that the pitching motion is about, and the 
pitching moments are referred to, the leading edge or apex. 


The position of the axis of pitching motion, and of the axis to 
which pitching-mioment coefficient is referred, is of importance. In 
this j?eport, these axes coincide in an axis normal to the root chord of 
the wing and x>as 3 ing through the wing leading edge or apex. However, 
it is often desired to transfer the pitching motion and mcment calcula- 
tion to other axes, and the formulas for such a transformation will be 
given here. Let the pitching motion refer to an axis lying a distance 
aco back of the leading edge or apex, and let the pitching-moment 
coefficients refer to an axis bcQ aft the leading edge or apex. The 
necessary transformation formulas are 


a ^ 

b ^ ~ ^ ~ 

where a prime on a quantity indicates that the pitching motion is about, 
and the pitching moments are measured about, an axis through the leading 
edge or apex of the wing. ^The quantity q is the dimensionless rate 
of pitching, equal to Cq 9 /Vq, The subscripts a and b mean, 
respectively, that the quantity in question refers to a pitching motion 
about an axis at \a distance ac'Q aft the leading edge or apex, or that 
the pitching moments are measured about an axis bcQ aft the leading 
edge or apex. 
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WING WITH SUPEESONIC EDGES 
INDICIA! LOADING FOE SINKING WING 
Analysis 


The seven regions .- The analytic expression for the Indiclal load- 
ing over the triangular wing has a different form In each of seven 
regions. These regions are determined hy the positions of the various 
wave fronts relative to the wing plan form (fig. 1 ). For t<0 the 
wing la motionless. Its leading edge lying along lines represented hy 
the dashed lines In figure 1 . At t =0 the wing starts suddenly to 
move, and for t> 0 , travels forward at a constant speed Vq, After a 
certain time t has elapsed, the wing has traveled to a new position, 
also shown In the figure. In this same Interval of time, pressure 
Impulses have traveled out In spherical waves from every point of the 
region which the wing has occupied. The trace on the wing of the 
sphere starting from the wing apex at t=0 forms the external boundary 
of region 7. The area outside this circle and within the traces of the 
cylindrical waves (the envelopes of the spherical waves) generated hy 
the leading edges at t =0 forms region 4 . Eeglon 5 Is formed hy the 
overlapping of these oyllndrlcal waves, and the solution for loading 
within It can he found hy a suitable superposition of the solutions for 
regions 3 and 4 . Eeglon 1 lies between the cylinder trace on the wing 
and the leading-edge position at time t; the loading In this region 
cannot he affected hy the manner In which the wing started Its motion 
since It lies outside the starting cylindrical waves. Bence, the 
loading In region 1 1s the same as that on a swept wing flying at a 
steady supersonic si)eed. The solution In region 2 can also he obtained 
from steady-state lifting-surface theory, but, whereas In region 1 the 
field Is two-dimensional (l.e.. Invariant with distance measured parallel 
to the leading edge), in region 2 the field Is conical. Eeglon 6 is 
formed hy the overlapping of regions 2 and 4 . Finally, region 3 Is that 
area completely unaffected hy waves from the wing edges. In the follow- 
ing subdivisions the analysis of each of the separate regions will he 
discussed. 

Eeglon 1 ; The loading In region 1 of figure 1 is equal to the 
loading on a two-dimensional flat plate moving at a constant velocity 
given hy the component of stream velocity normal to the leading edge 
of the triangular wing. Since this component is supersonic, the 
loading is' of the Ackeret type and is given hy 

f 

\ 3 “ 

But since = Vq cos A where A is the angle of sweep (see fig. 2 ), 
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cos^a 


Cjj = a sec A 
= Mq cos a 

=Jy^ cos^A-1 


and 



Finally, if ctn A = m 



Ijcon 

ypv - 1 


( 5 ) 


where 


3 =yMo2 _ 1. 


Eegion 2: The steady— state loading on a triangular wing with 

supersonic edges has teen given hy several authors (see, for conven— • 
ience, reference 5) so the expression for the loading in region 2 can 
he written immediately for the coordinate system shown in figure 2 as 



7T + arc sin 


3^my— (x+Mpt ) 
3[m(x+Mot)-y] 


arc sin 


3^my+(x+Mot) \ 
3 [M(x+Mot)+y] / 


(6) 


Region 3 ; Since region 3 is unaffected hy the edges of the wing 
the solution for the loading therein can he written as in reference 6 



(T) 


Region 4 ; The solution for loading in region 4 can he obtained 
from consideration of a two-dimensional wing starting from rest and 
moving with velocity normal to its leading edge. This problem 
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has been treated in reference 6 and the solution written there can he 
written for the right-hand side of figure 1 as 



_ ^n 

r 

MaXH+t . -s/Mn^-1 
arc cos — + — 7^ 

^ — + arc sin 1 

V Ini 

7T3n 

1 

L 

Xn+I^t Mn 

V 2 * y J 


where the notation, as defined by the sketch, is 


Xjj - X cos A — y sin A 
y^j = X sin A + y cos A 

and since 

1 


ctn A = m, sin A = 


\/ l+m^ 


cos A = 


m 


»J l+mP 


Then 


mx— y 

Xn = ■ 

V 1+n? 

x+my 

yn = 7= 




The equation for loading now becomes, in the coordinate system of 
figure 2, 



k<xm 

mMo(mx — |y| ) + t(l+H?) 
cos , - ... . + 


f __ _ _ 

4 It*/ 1 

■_rvX vy • 

Vl+n? (mx — |y 1 + mMot) 


U * arc Bln 

“Mo t-/IT^ / 


( 8 ) 


Region The solution for loading in region 5 can be obtained 
by superposition of the solutions for regions 3 and li-. If the solu- 
tions for the two sides of region 4 (obtained from equation 8) are 
added, the result gives a value of w^ of twice the required amount 
in region 5^ as well as undesirable pressures off the wing. However, 
subtraction from this sum of the solution for region 3 (equation 7) 
reduces the downwash w^ to the proper value, and also cancels the 
excess pressures. The resulting expression can be written 
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4am 


rt 


*[ 


mMo(BEC-y) + (l+iE?)t 

arc cos — — + 


‘/l+ii? (dbc— y + mMQt) 


mMo (mx+y) + (l+n?) t 
arc cos — ■■ + 

a/ 1+m? (mx + y + mMot) 


»J { . mx+y 

I arc sin — ^ — + arc sin — 

“Mq \ t 7 1+n^ t 


mx-y 


(9) 


Region 6 ; The loading in region 6 can also he calculated hy 
supei^osition. To find the loading in this case^ add the solutions 
for regions 2 and 4 (equations (6) and (8)) and subtract the solution 
for region 1 (equation ( 5 ))» There results 



4om 

■SitJ 02)1^. 




p2jny _ (x+Mot) 

sin — : — ^ • — arc sin 

3im(x+Mot) — y] 




(x+Mpt) ^ 
3 [m(x+Mot )+y] 


mMo(mx - | y| ) + t(l+m^) ^ 

(nix - I y| + mMot)<\/l+n? 



n 

2 


+ arc sin 


mx - |y| ) 1 
tA/l+nf ^ 


( 10 ) 


Region 7 ? The solution for the 
obtained by means of equation (3), 



loading in region 7 can be 

The analysis used in finding the 
solution in this region is 
not difficult but the algebra 
is rather involved. It is 
useful at this point to 
Introduce polar coordinates 
(see sketch) such that 

X— Xi = r cos 9 

y— 7 i = r sin 9 
dxidyi = r dr d0 (11) 

From equation (11), eqxxa— 
tion (3) can be written in 
the form 

a 
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The acoustic plan form for points in region 7 is "the region 
bounded by three curves as indicated in the sketch. The arc between 
01 and 02 ia determined by eliminating T between the eq.uatlQns 

r^ = (t - T)^ 

(the eq.uation for the inverse sound waves) and the equation for the 
left leading edge 

yi = -m(xi + MqT) 

The arc between ^ and 0 i is found by determining the acoustic 
intersection (eliminating T) of the right leading edge with the 
inverse sound wave; and, finally, the arc between 02 and 03 is 
given by the equation r=t. The equations of these arcs can be 
written in polar coordinates as 
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in(x+Mot )+y 
m cos 0+mMo+sin 0 


01 ^ 0 ^ 02 


r = t; 02 ^ 0 < 03 

r = m(xfMot)-y f q < 0 < 

m cos 0+mMp-sln 0 I 03<0 <2n 


Using these expressions, equation (12) reduces to 

(p = 1°^ m(x+Mot)+y ^ ^ 

J JO, cos 0+mMrt+sin 0 


01 

0 , 


t(30 + 


2n 


- f 
" '^2 

A m(x^■^ 

^ Ja COS 0 + 


m(x-i-Mot)-y 


jQ m ous y +mMQ— sin 0 


d0 


and taking the partial derivative with respect to t (to determine 
the loading according to equation (2)) one finds^ 


( ^') = 2a_m 

\ <lo^ ^ uu mMrt + m 


d0 


>0^ ioMq + m cos 0+sin 0 itHq 


2 g r 

Jn_ 


d0 + 


2a m r ^ 

« Jq 


d0 


(13) 


'0g mMQ + m cos 0-sin 0 

In evali^ting this equation^ the following integral is used! 
for — rt < 0 < jt 

f — = - ■ arc tan °>(Mo-l)tan (0/2) ±1 

idMq + m cos 9 ± sin 9 y i ' 

The limits di, 9 z, and 6q are all functions of t but in moving 
the .^rti^ derivative through the integral sign the terms involving 
o0i/ot, 002 /ot, and 803 /8t all cancel one another. 
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Since equation (lit) is valid only in the interval -n ^ 0 ^ n, care 
must be exercised in applying it because the angle 9 q my be greater 
than n (as in the preceding sketch). In case n< 0 , it is con- 
venient to introduce the angle 63 ’■ = 0g— 2 jt» The expression for 
^p/To cflu then be written in two forms, according as 0 is less 
than or greater than n: ® 


n + arc tan (9i/2)-l _ 

V P^a^-l 


for y > 0, 0g <Jt 

( ^ ) _ ^kxm 

^ ^7 « >/ p^n?-l L 


arc 


arc 


tan + arc tan ^(Mo-Dtan (9a/g)-H 


m(Mo-l)tan(03/2)-l 1 . 2a ^ ^ 


for y > 0, Jt < 0, 


. m(Mo— l)tan (0i/2)— 1 

1 =*• ) = — — — ' ■ arc tan — rr — — ■■ — 

V^lo^ n </ p2D^_i L 


^ -(Mo-l)tan (9i/3)tl ^ m(Mcr-l)tan (6^/2 )tl 

y p^n^-l yp^mS-l 


arc tan 


m(Mo— l)tan (03’/2)— 1 


* + 2") (15b) 


where 


0, = arc cos f Mpt ) ^(x^Mot )^-p V 


y2 + (x+Mot)‘= 


(0 < 01 < It) 


02 = arc cos 


m(y+mx) — (l+m^)t^ — (y+mx)^ 

(l+m?)t 


(0 < 02 < It) 
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= arc cos 


-m 


(y-mx) - ^/ — (y-nmc)^ 

(l+m?) t 


03»= 03 - 2n (if It < 03) 


The limitation on 0g can he given both an analytic and geometric 
interpretation. Thus, equation (15(a)) applies. for 0 < m(x+t )< y 
and equation (15(h)) applies for 0 < y < m(x+t). These regions are 
shown in the acconqianying sketch. Because of the. geometrical symnetry 
about the x axis, equations (15) suffice for the determination of 
loading throughout region 7» 



Discussion of Eesults 

Plots of the load distribution on the sinking triangular wing with 
supersonic edges are shown in figure and an isometric drawing of the 
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loading on the right panel appears 
in the sketch. The positions of 
the spanwise sections were chosen 
so that each of the regions 1 
throvigh 7 is represented. It is 
to be noted that the results for 
region 7 show no unusual charac- 
teristics, but fit in well with 
those for the adjoining regions. 

In general, the distribution is 
similar to the steady— state load- 
ing on a triangular wing. 

UroiCIAL LIFT AKD 
PITCHINO MOMENT 
Analysis 

Methods of solution .— If a 
detailed knowledge of the load 



distribution is not required, but onl;y the total values of lift and 
moment need be known, analyses much simpler than the one presented in 
the previous section can be employed. These methods, however, not only 
require the edges of the wing to be supersonic, but also require the 
trailing edge to be straight and normal to the free— stream direction. 
One such simple method has been presented in reference 4 , It involves 
the integration over the three-dimensional plan form of longitudinal 
strips, or elements, which carry the two-dimensional values of loading 
as a function of time. The results presented in the present section 
can be derived by this method as well as by the method to be developed 
next. 


Consider again equation (l) and integrate each term with respect 
to y between the limits minus and plus infinity.^ There results the 
equation 



If y = and y = yj.(x,z,t) are the equations of the Mach 

waves streaming back from the leading edgss on the left and right sides 
of the wing, respectively, then, since 9 is continuous across these 


2 — — ■ ■' 

'The basic idea for this solution was given by Prof. P. A. Lagerstrom 
in his lectures a t the California Institute of Technology. 
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waves tut and are not. 


^yr 




89 


8 r'\. ^yr ^yz 

8x® J 8x 8x J 3x^ 

y, . 7i 


cLy 


where Uy. and U2 are the values of u on the interior faces of the 
right and' left Mach waves, respectively, and 




£1 

8y^ 


cLy = - V, 


Values of the terms involving w and 9 ^ are similar to those involv- 
ing u so that finally, if 


then 


$ = 


9 dy 


y-, 


(16) 


2^2 
8 ~^ 8 <ji _ 8 $ 

8x^ ^ 8z^ -8t^ 



8 yr 

Uy. + Vy, - 

^yr 

\ 8 x 

8 z 

St 

r' ^yz 

8 yj 

^yz 

V 8 x 

^ 8 z ‘ 

St 


^t 


-vr)- 


■Pt 


The terms enclosed within the brackets in the last equation combine 
so tliat each bracketed quantity is zero. For the case of interest here, 
this is not difficult to show. Consider, for example, the right wedge. 
Then, since the equation of the wedge is 


yr 


= — z »/ — 1 + 


mx + mMot 


and the value of the potential is the steady— state two-dimensional value 
given by the expression , 




m(x+Mot )— y— z p^n^— 1 
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the term 





"becomes 

-Wy ( m — — + J 32bi?-1 

g 

and this is identically zero. 

Finally, therefore, eqixation (l) has "been reduced in terms of eqxia— 
tion (16) to 

ht - ®xx - «zz = 0 (IT) 



Consider next the "boundary conditions for a triangular wing with 


supersonic edges. In the plane z= 0 , 8$/8z "becomes 


w^dy where 


w^ is the vertical induced velocity in the plane of the wing. Since 
the edges are supersonic, however, w^j must he zero off the plan form 
of the wing and the integration need extend only over the plan form 
itself. Hence 



= -^YqCoh (x+Mot) 

z=0 


(18) 


for the flat triangular wing at constant angle of attack and 

^1^') =.-29m (x+Mot)^ 

Adz/z=o 


(19) 


^It is not necessary to perform a direct calculation in order to prove 
the above result for arbitrary plan forms. The bracketed terms 
represent the directional derivative of the velocity potential taken 
along the so-called "conormal" of the foremost distiirbance surface. 
Since cp is constant on the surface, and since the conormal lies 
along the surface, the bracketed terms are zero. 
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for a flat triangular wing with a constant rate of pitch 9 about its 
apex (positive 5 produces a downward motion of the trailing edge . 
relative to the leading edge). 



These boundaiy conditions are exactly like those studied in steady- 
state supersonic wing theory. In fact, the lifting-svirface analog 
(shown on the sketch) is a wing tip of speci fied camber in a supersonic 
free stream having a Mach nunber equal to 2. The solution for the 
potential in the plane of the wing for this problem can therefore be 
Witten immediately as 


( f ) 


z=0 


1 

jt 



8z /z=o 


*/ (t-ti)^ - (x-xi)^ 


dti dxi 


where a is the portion of the area on the shaded surface in the sketch 
lying ahead of the forecone traces given by the equation (t— ti)^ = 
(x-xi)^. Using eqmtion (2) for the loading coefficient, . and introducing 
the following notation for the average spanwise loading 

^ J_g (sinking wing) (20a) 

Pi = ? f ^-dy (pitching wing) (20b) 

23^ 0 ^-s <lo 
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where 


8 = m(x+MQt) 


it is found that 

P = hE A. f r dxidti 

° jt3 8t ^ J y (t-tx)2 - (x-xi)2 

^ " rts2 at ^ J ^(t_tx)2 - (x-xx)2 


( 21 ) 


( 22 ) 


(23) 


Ving with constant angle of attack (sinking wing),— The solution 

for the average load on a flat, supersonic-edged, triangular wing 
starting from rest at t=0 and flying at a constant speed and angle of 
attack is given hy equation (22). With the transformations 


X— Xx = Xe 


‘b— t X = \ j 2. 


this hecoiDes 


Po 


Urn a 

Its at 



dtgdxE 


8 »» 4 . 

m - xg - Mptg 

V ts^ - X2^ 


This integral can he evaluated and gives for x<t (region A) 



(2i*a) 


for — t < X < t (region B) 


'■o = ? 




^ + 


x+Mjjt 

and for x < — t (region C) 


arc cos 


+ Mo 

\ V 3 


arc cos 


t+Mpx 

x+Mot 


(24b) 




(24c) 


where the regions are shown in the preceding sketch. Equations (24) can 
also be obtained by integrating the equations for the loading given in 



22 


NACA TN 2387 


the preceding section. These integrations were carried out (in some 
regions numerically) and the results were found to agree with those of 
the present analysis. 

It is now possible to write the equation for the indicial lift and 
pitching moment for a sinking wing: 

* 

.Co-Mot 

Ct = — ^ / m(x+Mot) P dx (25) 

Cm ' = - r — ^- 7 T / m(x+M t)^ Po dx (26) 

a . S Co Mo 

where S is the wing area (equal to mCo^) and the prime indicates that 
the pitching .moment is measured about the apex, the positive moment 
being one which causes the trailing edge to sink relative to the apex. 

Combining equations (24) and (25) one finds for the first interval 
shown in the sketch 


Interval no. 
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‘^Mo 


(x+Mqo ) dx + 


1_ x+Mq^ 


+ arc cos 


:-f) 


Mr, t+MoX 

arc cos — , 
p . x+Mot J 


dx + 


Co-Mot 


4(x+Mot) dx 


This eqmtlon integrates to give, if to = t/co, 

for 0 < to < — — (first interval) 

- ° - Mo+1 


Or = iL (1 + 1 to^ 

Mo \ 2 ° 

Similarly for ;■ ^ to < ■ (second interval) 


Mo+1 


Mo-1 


(27a) 


/ 


1 + 2 


=377- (1 + r to Jarc cos 


-'a . jtMo 2 y - t 


Mpto-l 


^ arc cos (Mo~to3^) + 


2 ^2^ - (1 - Moto)‘ 


(27b) 


and for 


1 

Mo-1 


^ to 


(third interval) 


Cl 


a 


U 

3 


(27c) 


In the same manner the values for Cj]^* in the various intervals can be 
determined by combining equations (24; and (26). There results 
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for 0 < < 


Mo+1 


(first interval) 




3 Mo 

for ■ ; ^ to < — i — (second interval) 

Mq+1 M -1 


-) 


(28a) 


% jtMo 


i (8-Moto-Mo^to2-2to®) ytoMl-Moto)^ + 


^ ( 2 +Moto®) arc cos ^ arc cos (Mq— P^ to) 

t© 


(28b) 


and for -r^- ^ to (third interval) 
Mq-1 

C • = - 8- 




3 P 


(28c) 


Numerical results for a Mach number of 2 will be presented in the dis- 
cussion section. 

Ning with linear angle-of-attack variation (pitching wing) .- The 

solution for the average load on a flat, super sonic-edge d, triangular 
wing flying at a constant speed and pitching at a uniform rate § 
about its apex is given by equation (23). With the transformation 


X-Xi = Xs 


t— t 1 = tg 


this eqiiation becomes 


P-. = 


kvF 

JtS^ 


ill 


- - X2 - Mota 




dta dx2 
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and the evaluation of this gives (for the regions defined for the 
sinking wing) 

for X < t (region A) 


Pi = 



^Y1 

x+Mot/ J 


(29a) 


for — t < X < t (region B) 



’ 1 t^+jx+Mptf 
_ 2 (z+Mot)^ 


arc cos 



+ 


Ma 


arc cos 


1 3k+2M^t 
x+Mot 2 (x+Mot)^ 


(29h) 


and for x < — t (region C) 



(29c) 


The equations for lift and pitching moment on a pitching wing can 
he obtained from the equations 


Cl’ 


(0 


= Cl 


* = — / m(x+Mr)t)^ Pi dx 

SCO Mo 


(30) 


Cm* 


- Cm^’ - 


^ f 


°o“^o'*" 


SCo*^Mo 


^t 


m(x+Mot) Pi dx (31) 


where the primes indicate the wing is pitching about and the monents are 
measured about the leading edge. 

A combination of equations (29) and (30) gives for the lift coeffi- 
cient 


for 


0 < to 


<_J_ 

Mo+1 


(first interval) 
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(32a) 


for . ■ ^ t ^ . (second interval) 


Mo+1 


Mo-1 


Ct ’ = 


8 \ f 1 .. z Iw 3 1 ^ 

s; L U ° ■ 3"°*® * 3 / 


Motn-l 

arc cos ^ ^ + 
to 


33 


arc cos 


(Mo-3^to) + (^ I t^2 ^ 1 Mo^to^- ^ Moto + 
v^to^-(l-Moto)^ ] 


(32t) 


and for — i— < to (third interval) 

Mq -1 


Cl^' = 8/33 


(32c) 


Similarly a comhination of equations (29) and (3I) yields, for the 
pitching moment about the apex, the results 


for 0 < t_ < 


Mo+1 


(first interval) 




Cmq' ^ ~ ^ (1+^*^^) 


(33a) 


for y < to ^ Fi^r interval) 


Mo+1 


C«, ' = 
^Diq 


«Mo 


2(l+to^) - ^ V (1+H)^) 


Moto-1 

arc cos — r + 


^ arc cos (Mo~3^to) + 


12 


42-22Moto+ (2Mo^+3 )to^+ (2Mo^+13 )Moto° 


/'Ao^~(l-Mo'to)^ ^ (33t) 
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and for 



(third interval) 



(33c) 


Numerical, results for a Mach muuiber of 2 are presented in the 
next section. 


Discussion of Results 


Figure 4 shows the variation of location of 

^c.p. with the number of chord lengths traveled ( Tq = /cq = 

Mot/co = MqIq) for a free— stream Mach number Mq eq.ua! to 2. The 
values of C]^ and Cn^' are in agreement with those given in refer- 
ences 1 and 2. For the purpose of comparison, similar curves are shown 
for a two-dimensional wing having a chord Cq equal to the root chord 
of the triangular wing, and for the triangular wing in reversed flow. 
The material necessary for the calculation of the latter curves was 
presented in reference 6, Several interesting conclusions can be drawn 
from these results. 

First, notice that the total indicial lift on the triangular sink- 
ing wing is the same at every Instant as that on the sane wing in 
reversed flow (both wings of course having started with the same velo- 
city at the same time), and that the value of this lift is the same as 
the total indicial lift on the two-dimensioneil wing only at the begin- 
ning of the motion and again when the steady state has been attained 
(fig. 4(a)). Such a result for reversed flow is not true for the 
pitching moment (fig. 4(b)) and center of pressure, and it can be shown 
that it is true for the total lift only when the wing is a flat lifting 
surface with supersonic edges^ 

Second, notice that, since all of the characteristics for the 
triangular wings are Independent of the angle of sweep, they are valid 
for any unyawed triangular wing flying at a Mach number equal to 2 
and having supersonic leading edges. 

Third, it is apparent from figure 4(a) that the transition of the 
total indicial lift from its initial to its final value is less abrupt 
than that transition for the two-dimensional wing, and finally the 
movement of the center of pressure on the sinking wing of all types is 
seen from figure 4(c) to be small. 


R. T. Jones has shown, in an unpublished work, that the build-up of 
total indicial drag on symmetrical nonlifting wings is the same for 
all types of plan forms in forward and reversed flow. 
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Finally, it can te seen from equations (27a) and (27c) that 
the initial and final values of depend on l/M© and I/3, respec— 

tively. As the Mach number Mq is increased, therefore, the variation 
dies out since P and Mq become neeirly equal. The same remark 
applies to all the other coefficients. 


Figure 5 presents the values of , and Xc.p./°o 

wings pitching about the foremost extremity of their plein form (i.e., 
leading edge or apex). Again the results are presented in terms of 
chord lengths traveled ' Tq for a free— stream Mach number equal to 2. 


For these wings it is apparent that the reverse— flow theorem does 
not apply even to total lift. The results for the triangular wing are 
still independent of the angle of sweep, however, and the movement of 
the center of pressure is again slight. 


WING 'WITH SUBSONIC EDGES 
INDICIAL LOADING FOR SINKING WING 
Analysis 


The six regions .— As in the study of supersonic-edged triangular 
wings, there are also in the case of triangular wings with subsonic 
leading edges various regions in which the analytical form of the load- 
ing equation is different. Figure 6 shows the regions into which the 
subsonlc-edged triangular wing can be most conveniently divided. Most 
of these regions have counterparts on the supersonic-edged wing shown 
in figure 1. 

To begin with, region 6 lies within the spherical wave which 
started at t=0 from the wing apex. Region 1 is within the cylindrical 
wave which was started at t=0 by one edge of the wing, but outside 
the wave started by the other edge. Region 4 is the area formed by the 
overlapping of the two cylindrical waves from the opposite edges, but 
outside the region influenced by the reflection of one of these waves 
on the opposite edge (secondary wave fronts shown in the figure). 

Region 5 is the area between regions 4 and 6 where the flow is influ- 
enced by secondary (and higher— order ) wave reflections. Finally,' 
regions 2 and 3 are similar to regions 2 and 3 in the supersonic-edged 
case; region 2 being that uninfluenced by the starting phenomena and 
therefore having a loading already at its steady— state value, nnH 
region 3 being that which is unaffected by the disturbances emanating 
from the edges. 
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Region 1 ; The solution for the load dlstrihution in region 1 

is the same as that for a two-dimensional wing starting suddenly from 
rest and moving with a steady subsonic velocity Vn. A solution to 
the latter problem for the initial part of the motion is presented in 
reference 7» In terms of the normal components of velocity and dis- 
tance, therefore, the loading coefficient for the ri^t-hand side of 
figure 6 can be written inmediately: 

(jh- 

V'ln /i ’^Vn^n \l+% v ^^t+xn V t-Xn / 


The equations which relate the normal components to those in the 
free— stream direction have already been given in the section on 
region 4 of the supersonic-edged wing. Use of these relations leads 
to the following expression for loading in region 1 (in the coordinate 
system of fig, 7) 

f ^ + |y| — mx ^ 

V qo ^ \ iiiMq+ ^i+np tj m Mot — | y | + mx 


arc tan 


mMot — lyl + mx 
t Vl+n^ + |y| — mx 


m 


Region 2 ; The loading on region 2 , being the steady— state 

loading on a triangular wing with subsonic edges, is well known. The 
solution for region 2 of figure 7 is therefore givenvby (see, for 
convenience, reference 5) 

( 


= 

^0/2 


4anP (x+MQt ) 

''/mP (x+MQt )^— y^ E 


( 35 ) 


where E is the complete elliptic integral of the second kind with 
modulus 

Region 3; The loading in region 3 follows from reference 6 and 


is 
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Beglon 4 ; The loading in region 4 of figure 7 is calculated hy 
superposition. Just as the solution for legion 5 of the wing with 
supersonic edges was obtained. The solution in region 4 is the sum 
of the solutions for the right and, left halves of region 1 , minus the 
result for region 3« Thus 


<lo/4 


8 a 


mMo 


L H1M0+ <v/l+D^ 




t «/ 1 +iif -y-mx 
mMot+y+mx: 


tV l+#+y-mx 'N 
mMot— y+ncc J 


+ arc tan 


A 


mMQt+y+mx; 


l+n?-y-mx 


arc tan 


mMot — y + mx 


t i/l+m? + y — 


mx 


-1] 


(37) 



apparent that the principal difference 
the behavior around the leading edges; 
supersonic edge, whereas it becomes infinite at the subsonic edge 
view of the known steady— state results this difference was to be 
expected. Elsewhere the loadings are q^uite similar. 


Regions ^ and 6 ; In these 
regions the exact solution for 
the loading has not been deter- 
mined. As was shown in refer- 
ence 4, such solutions would 
require the solution of a three- 
dimensional elliptic— type partial 
differential equation. In this 
report a later section will con- 
tain an approximate solution for 
these regions. 


Discussion of Restilts 


An isometric drawing of the 
load distribution, for the 
regions in which it is known, is 
shown in the sketch. Comparing 
the results for the loading on 
this wing to the one with super- 
sonic edges (fig. 3 ), it is 
in the two distributions is in 
the loading being finite at the 

In 
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The results presented in equations (3^) through ( 37 ) will next be 
examined in a different light. Choose a given spanwise section on the 
wing and watch this section as time progresses from t=0. This amounts 
to fixing the axis on the body and can be accomplished simply by using 
the quantity s introduced in equation ( 21 ), 

s = m(x + Mot) 

It is clear that s is the semispan of a given spanwise section, and 
that if equations ( 3 ^) throu^ (37) are written in terms of s, y, and 
t, for a fixed s they represent the variation of loading on a given 
section as time progresses. 

If the notation is further simplified by introducing the parameter 
Pq where 


P = 7 

® mMo+ »f L+n? 


( 38 ) 


equations (3l)-) through ( 37 ) can be written in the following way: 


/IKIE 4 IM * aro tan / , ° - 1^1 \ 

1 \ S - |y| t/Pg - 8 + |y| / 


T— ) = 


4ams 


E V s 2 — y 2 


Ap \ _ 4a 

Ho A Mo 


\lo'3 


(39) 

(40) 

(41) 




- y - s 


arc tan /- - - + arc tan / ? '*' ■ I ... - i j (42) 

V t/Pe + y - s V t/Pg - y - s 2 / 


The load distribution anross any section is given by eqmtions (39), 
(4l), and (42) from the time t/Pg = 0 to (t/Pg)i, where the term 
(t/Pg)^ is equal to 2s or s/m(Mo+l)Pg, whichever is smaller. 
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(At t/Pg = 2s the spanwise section has just reached the secondary 
waves shown in figure 6, and at t/Pg = s/m3g(Mo+l) the spanwise section 
has Just reached the spherical wave which started from the apex. ) From 
(■t/3e)i ~ loading has not been determined, 

and from t/Pg = (t/Pg)^ to t= » 
the loading is the steady— state 
value given by equation (40). 

The sketch shows this initial and 
final load variation plotted as a 
function of the parameter t/3g. 

At the beginning of the motion 
the loading is constant across 
the span, but this type of dis- 
tribution is quickly modified and 
the shape of the curve tends 
toward the "inverted elliptic” 
loading given by equation (ij-O) 
and shown in the sketch as the 
distribution at t/pg = (t/Pg)^, 

In fact, when the span has 
traveled a distance such that 
t/Pg=2s, the expression for the 
loading given by eqxiation {k2) 
becomes 





16 mPgS 

2s It */ s^— y^ 


(h3) 


which differs from the value given by equation (^l-O) only by a constant 
of proportionality. Both before and after the time t/Pg=2s the shape 
of the loading curve varies from the simple type represented by equa- 
tion ( 43 ), but the trend, and to a certain extent the rapidity of the 
trend, is clearly established. 

The average spanwise loading Pq, introduced by equation (20a), can 
now be determined for certain regions. Hence, if the notation 


T = t/s 



( 44 ) 


is adopted, there results for the early part of the motion, that is, for 

0 < t/Pq < ( t/Pq)i 

Pg = 2 (^2 - ^ ^ + 4mPgMo ( T/Pg) (U5) 
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Eq.mtion (45) vas derived by integrating equations (39), (4l), and (42). 

For values of , v/Pg > 1 /Pq]u(Mo— 1 ) equation (40) is valid. 

Hence 



P 


o 


2wmMo 

E 


(46) 


The sketch indicates the 
magnitude of this average load 
for both large and small values 
of ''■/3e» Notice that for small 
values of x/p^ it is suffi- 
cient for the establishment of 
the curve to specify the parameter 
mMoPg, but for large values an 
additional parameter must 'be given 
(such as Mq in the sketch). 
Notice, further, that in spite of 
the large variation in the distri- 
bution of the loading, as shown in 
the previous sketch, the average 
value Pq varies linearly through- 
out the Intervals considered. This 
result is similar to the one 
obtained for triangular wings with 
supersonic edges and is given in 
equations (24). 



INDICIAL LOADING ON VERY SLENDER TRIANGULAR WINGS 


AneJ-ysis 


In the first section of this report entitled "The Problem and the 
Nature of its Solution," it was pointed out that if the wing is slender 
(i.e., has a small ratio of spein to chordwise length) the basic partial 
differential equation (l) can be approximated by the equation which was 
-previously introduced as equation (4), thus; 


9tt -9yy -CPzz =0 

The boundary conditions appropriate to this problem will now be examined 
in some detail. 
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Just as in the previous sections of this report, consider a triangu- 
lar wing which is at rest for t < 0, starts suddenly to move at a for- 
ward velocity equal to Vq at t=0, and continues at this same velocity 
for t > 0. It should he emphasized that in this case, Vq may he either 
subsonic or supersonic* A section in the spanwise direction, as for 

instance section AA in the sketch, 
projects into the yt plane as a 
narrow rectangular strip along the t 
axis. Since equation (4) has heen 
derived on the assumption that the 
velocity gradients in the y, z, and 
t directions are independent of the 
gradient in the x direction, the 
hoimdary conditions along the strip 
shown in the sketch are independent 
of those on other strips projected 
from spanwise sections along the wing. 

I ‘ Hence, the problem is to find a solu- 

tion to equation (4) which will make 
cp^ constant over the strip and at the 
same time will satisfy the other con- 
ditions listed under equation (l). 

In the lifting— surface analog this 
corresponds to the problem of finding 
the velocity potential over a flat 
rectangular wing of low aspect ratio 
situated in a free stream moving at a 
Mach number equal to Solutions 

to the latter problem can be obtained 
by various techniques, and so the pro- 
cedure will be first, to find the potential for the steady-state, flat, 
rectangular wing, and then, by analogy, to convert this to the solution 
for the very slender triangular wing in unsteady motion. 

The steady-state, lifting-surface problem .- Lifting-surface solu- 
tions for the loading on a rectangular wing traveling at supersonic 
speeds have been developed for regions 1, 2, and 3 of figure 8 (by 
Busemann and others), and by means of these solutions the load distri 
button on a spanwise section of the triangular wing can be determined 
bo a time necessary for sound to travel that span length. For t > 2s, 
however, the solution becomes considerably complicated by the increasing 
number of reflections from the, edges. Heference 8 gives solutions for 
the loading on a rectangular wing in region 4 and indicates methods for 
extending the solution to regions farther along the wing. Already in 
region 4, however, the expression is cumbersome and in higher-numbered 
regions the expressions become difficult to manipulate. These methods, 
therefore, will be discarded in favor of a more approximate but simpler 
analysis. 

If X is the distance along the chord, y the distance along the 
span, and s the semispan, then the solution for regions 1, 2, and 3 of 
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figure 8 can be written {far conTenience, see reference 5)t 


Region 1 



Region 2 


Region 3 


^o 


8a , 

— arc tan 

jt 




^ _ 8a 


(' 


arc tan 



3+y 

x-s-y 


+ arc tan 



g-y 

x-s+y 



(i^7a 


(47b 


(47c 


As X increases (i*e.j for higher-nvuabered regions in figure 8) 
it is reasonable to assume that the spanwise variation of loading is 
relatively unimportant — except that it be "smooth" and fall to zero at 
the side edges — and the chordwise variation of loading is domiriant. 
Assume, therefore, that the loading is given by the relation 
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Spanwise^ this has the variation shown in the preceding sketch; chord- 
wise it is as yet arbitrary. To fix the chordwise distribution the 
value of f(x/s) will be determined so that the vertical induced velo- 
city along the center line is constant and equal to — Voa. 

The solution to this somewhat artificial problem approaches the 
exact solution to the steady— state lifting— surface problem for a flat 
rectangular wing along sections far behind the leading edge; closer to 
the leading edge it only approximates the exact solution; and, of coiirse 
in the vicinity of the leading edge it will be least representative. 

But, on the other hand, the exact solution is known in the vicinity of 
the leading edge and it turns out that the solution of the problem posed 
above forms a reasonable continuation over the remainder of the wing. 

The velocity potential for the problem which has been set can be 
readily expressed in terms of an integration of elementary horseshoe 
vortices over the plan form. Since the Mach number equals ,^2, then 
according to reference 9, 

cp = r f (^/Iq) dxidyi. 

A (yi^+z2) J (x-xi)^ -yi^-z^ 

where A is the area on the wing within the forecone from the point 
P(x,y,z), at which cp is to.be determined (the shaded area in the 
sketch). . . 


The simplification of the last expression is given in reference 9. 
The result is the integral equation 


1 = f(T)) + I f f(.Tli) G )dq 

do ■ . 

where q=x/s eind G is given by 

E, 


G(q— qi) - 


Ea-d-ka^ )K^ 
ka 


1 ^ q-^, 


q-q^ < 1 


(h9) 


ki = ^ ; ka = T]--q^ 


The modulus of is ■ ki and the modulus of and Eg is ka. 
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The solution of equation (49) for f(ri) is not difficult when 
numerical methods are used. For intervals of equal to 0.2, the 

result is given in tabular form in appendix A, and also in the sketch. 



This result can be in^jroved in the interval 0 < tj < 2 by means of the 
correct solution in that interval given by equations (4-7). Comparison 
of the spanwise average of the loading given by equation (48) with that 
which can be derived from equations (4-7) gives an equivalent f(T|) in 
the interval which, when used in equation (48), will give the correct 
value of the average span loading. The sketch also shows a cui*ve for 
this equivalent f(T)) which starts at 4/it and falls linearly to zero 
at q=2. 

By using the sketch, or tte results listed in appendix A, the 
loading over a low-aspect— ratio rectangular wing flying at a Mach number 
equal to can be estimated. Of particular interest is the damped* 

oscillatory nature of the load, falling to zero at one span length 
behind the leading edge and taking alternately negative and positive 
values beyond this point.® A somewhat different approach to this 
problem (reference 11) has recently led to a solution very like the one 
given here. 

The unsteady analog, sinking wing .- The first step in deriving the 
imsteady-flow results for the sinking wing from the steady solution is 
to replace x with t. In eqiJation (49) this corresponds to replacing 
T) with T where t is equal to t/s (equation (44). The second step 
is to rederive the expression for loading coefficient since in the 
time— varying problem it is expressed in a somewhat different manner than 
in the steady— state analog. In the unsteady case, as the triangular 
wing moves through a fixed reference plane the local span intersecting 

■g— ^ — 

A discussion of this aspect of these results is contained in refer- 
ence 10. 
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this plane grows as a function of time and equation ( 2 ), which repre- 
sents the partial derivative with respect to time with x fixed, must 
he expanded to the form 


^o VqMo L 


8 acp 

2 



8 acp 

88 ^ 

_ 8 t J 

x~ 

\L ^ J 

s 

8 s 

t 


where 


8A9 

and 

8 Acp 

_ 8 t _ 

s 

8 s _ 


indicate derivatives taken at constant s 


and t, respectively. Since s is equal to m(x +MQt), 8 s / 8 t equals 
iuMq, and there results 


4o VqMo 


( 


~ 8 acp ' 

8 t 


+ mMo 


8A(D 

8 s 


D 


( 50 ) 

8A9 ] 


In the steady— state problems an analog to the term involving , 

-8s Jt 

is missing, and the loading coefficient is given entirely by an opera- 

r8A9" 


tion equivalent to 


8t 


It is necessary, therefore, to operate 


-“s 


further on the solution given for the loading in the steady-state 
problem to obtain the solution for the loading in the unsteady problem. 
But 


dti 


8Acp 


8A(p 

_8s J 

t '-'0 , 

- 8t _ 


so that if the notation 


~8A(P 

8t 


= Ia 


is adopted (where ( ^ J represents the loading in the analogous 

Vo/s 

steady— state problem), then the expression for the unsteady loading can 


be given in terms of 


Ap _ 1 
Q.O ^ 



by the equation 


[(t). (S 


*0 



NACA TN 2387 


39 


By the application of equation (5I) to equations (47), the loading 
for the various regions of the unsteady wing can he found. For region 2 
there results 


^ = 
<1« 


o >- 


8a 

— arc tan 
n 




■t-e + !y ! ^ da 



3-|y| 


4a dti + 



^ arc tan /IBIE dt A 
It V ti-s+ |y| y 


which becomes 


_ 8 ti 

qo ~ 


mMo 


/t— s- 
V s— 


t-e+ lyl 


+ arc tan 


/ 


a- lyl 

t-a+ |y[ _ 


The loading coefficient can he similarly derived in the other regions so 
that finally 


Region 1 


Ap 4a 
0,0 Hd 


(52a) 


Region 2 


_ 8a 
q© ■ itMo 


mMo 


' t-6+ 
S- I 


+ arc tan 


/ -"- | y | ^ 

y t-s+ |y I ) 


(52h) 


Region 3 


<lo 


8a / 


mMo 


^y ^t— s- 


s+y 

•y 


+ mMQ 


A-e-y 

V a+y 


+ 


arc tan 




a-y 

t^+y 


+ arc tan 



3+y 

t-s-y 



^ (52c) 
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For the interval t >2s eqiiation (48) must te considered. By means 
of equation (51 )> the expression for the loading coefficient can be 
■written 

^ (b)A& 


which becomes 


I0 Mq L 


(l-mMox) /l-(^) f ( t) + 
s 






r f(Ti)dTi'i (53) 

'^o J 


where f(r) is the solution to equation (49). Notice that for large t 
(when the loading has reached its steady state), f(r) is zero and 

/ f(Ti)dTi is unity. (See appendix A.) Hence, the loading is given 

Jo 

by the equation 


^ 4ams 


q.0 y s2-y2 

which is the steady^tate value for a slender triangular wing (equa- 
tion (4o) when E = l). 

It is now possible to derive the average span loading Pq as 
defined by equation (20a), thus. 


2sa J 9.0 


Placing equations (52) and (53) in this expression, it is found that 
for 0 <T< 2 

Po = 2(2-t) + 4mMoT (54a) 

and for x > 2 


Pq = «(l-mMoX)f(r)+2 mM^n f(xi) dxi 


(54b) 


Since the values of Pq given by equations (54a) and ,(5^t) were 
derived using different methods, their magnitudes at t = 2 are not 
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eq\Jal. The final curve for Po vsnat he constructed by fairing the solu- 
tion for T ^ 2 into that for x >2. The accompanying sketch shows 
these results together with the final curve chosen (solid line). 



The Tonsteady analogy pitching wing.— When the wing is pitching at a 
steady rate about its apex^ the equation for the vertical induced veloc- 
ity CHI the plan form is 

Wu = — (x + Mot)0 


so that the a in the steady— state equations (47) and (48) must be 
replaced by 0s/mVo» Since the loading coefficient is still given by 
equation (51), there results for the conversion of equation (48) the 
expression 


Ap 40 
^ - 


*/ 3^— f(''’) + mMo 


as 


^t/a 

3 3^-y^ I f(Ti) dXj 


axid this can be reduced to the form 
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^ ^ 
mMoVo 


(1-^dMot) */ s2-y2 f(T) + 


lj£ /2s2-yg 

\a/ s^— y^/ v7q 


f(Ti) dTi 


(55) 


As in the discussion of equation (53), it can he seen that equation (55) 
hecomes for the steady state (t large) 

^ ^ /2s^-y^ \ 

and this can he shown to agree with the steady— state slender-wing results 
given in reference 12, 

It is now possible to derive the average span loading as 

defined hy equation (20h) 

2b^9 ^-b *1o 


Using equation (55), one finds for t >2 


= n(l-^nMoT) f(r) + 3mMon / f(Ti) dr^ 


(56a) 


and a similar analysis based on equation (47) yields 
f or 0 ^ T < 2 

Pi = -^ (2-T+4mMoT- i mMoT2) , (56b) 

As in the case for Pq, the two equations for Pj_ do not Join at t=2 
qr|(i the final curve must be constructed by fairing the solution for t ^2 
into that for t > 2. , 
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Discussion of Eesults 


It is now possible to 
assess the accuracy of the 
solution for very slender 
wings in the Interval 0 ^2 

by comparing equations (52) 
and (54a) with equations (4l), 

(39), (42), and (45), the 
exact solutions for this region 
derived in the preceding sec- 
tion, It is apparent that the 
approximate solution differs 
from the exact only by a 
stretching factor in the t 
direction. Hence, if t . is 
replaced by and m 

(note m is proportional to 
y/t) by mPg, where pg is 
given by equation (38), then 
eqviatlons (52a), (b), and (c) 
are identical with equa- 
tions (4l), (39), and (42), respectively, and, of course, equation (54a) 
corresponds to equation (45), 

This rather remarkable result can be enlarged upon from another 
viewpoint. Suppose that in the steady-state analog problem the wing had 
been, flying at some Mach number other than say Mg, The solution 

to such a new problem could be obtained from the old one merely by 
applying the Prandtl-Glauert correction, that is, by stretching all dis- 
tances in the x direction by the factor l/Pg where Pg^= 

Such a procedure would convert, for example, equation (5^) to the form 

Po3g , = 2 (2- T /Pg ) + 4mPgMo t/Pq 

Finally, if Pq is adjusted so that Pq=^ at t =0, there results 

Pq = 2(2- T/Pg) + 4mMoT 

which is exactly the answer given by equation (45), It is possible to 
simplify the statement of this procedure by simply remarking: The exact 
resiilts for Ap/qg or Pq in the interval 0 < t/p^ < 2 can be 
obtained from the approximate results for a very slender wing by making 
an effective Mach number correction to the right-hand side of eqxia- 
tions (52) or (54a), respectively. 
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Actual trace of 


It is interesting to pursue this concept even further. Consider a 
spanwise section of a triangular wing as time increases from the starting 

impulse . The primary 
wave fronts emanating 
from either side pass 
across the section 
forming the Mach lines 
in the steady— state 
rectangular-^ing 
analogy. For very 
slender wings these' 
lines make a 45° angle 
with the trace of the 
side edge and are used 
to divide the plan 
form into regions as 
in the sketch. Now 
find the actual posi- 
tion of these primary 
wave fronts as they 
, form a trace on the 
section in the yr 
plane . A strai^t— 
forward calculation 
shows that these lines 
actually make an angle 
equal to arc tan l/Pe 
with the trace of the 
side edges. Hence the 
effective Mach number 
which is used to- 
correct the slender— 
wing results in the 

interval 0 ^ t ^ 2p@ is that which makes the Mach lines of the steady- 
state analogy coincide with the act\oal trace of the primary wave fronts. 



IHDICIAL LIFT AND PITCHiroJ MOMENT ON 
VERY SLENDER TRIANGULAR WINGS 

Analysis 


The lift coefficient for the sinking wing is given in the notation 
introduced in equation (21a) by the equation 


Cl = 


g 

WCr 


r 


2s 

Mo 


Podx 
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where Pq has been determined in the last section as a function of 
T=t/s. Consider the situation at a certain fixed time and let the x 
coordinate in the above formula be fixed in the wing. Then set 



( 57 ) 


and as before 


T = = mMo (58) 

^o 

where is the number of wing— chord lengths traveled. In this way 

the equation for lift coefficient becomes 

( 59 ) 

and by a similar analysis the pitching-moment coefficient taken about 
the apex can be written 


The equation for the lift and pitching-moment coefficients (where 
again the pitching moment is taken about the apex) on a pitching wing 
are 

of' 

and 


Discussion of Besults 


The values of Pq and P^ were taken from curves similar to the 
sketches in the last section (using the faired curves in the vicinity of 
T =2) and the results for the indicia! lift and pitching moment in terms 
of Tq^ the number of chord lengths traveled, are shown in figures (9) 
and (10) for a value of mMo equal to I/8, The results are all 
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q,ual 1 tat ively alike; in each case the curve falls from its high initial 
value at Tq= 0 to a minimum at about Tq=i/ 3 and then recovers and 

practically attains its 
asymptotic value at Tq= 1, 
This behavior is similar 
to that for a two- 
dimensional wing flying 
at subsonic speeds. For 
the piirpose of such a 
comparison the first part 
of the curve of indie ial 
lift coefficient versus 
chord lengths traveled on 
a two-dimensional wing . 
flying at a Mach number 
eq.ual to 0,8 (see refer- 
ence 7) is shown in the 
sketch. Notice again the 
rapid fall from the 
initial peak (Cl^=4/Mq) 

0 4 ^ 8 to a. minimum around 

^ To= 1/2, and then a smooth 

recovery to the asympto- 
tic value (Ct =2n/p), 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif., Jfar. 26, 1951. 
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APPENDIX A 

SOLUTION OF AN INTEGRAL EQUATION 

For convenience in applying equations ( 54 b) and (56a) of the text, 
a table of values of the function f(ri), obtained by nurerical solution 
of eqiaa-tion ( 49 ), is given here. The valioes of the integral, 

4 f(tli)dTji y are also listed. 
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wing with supersonic leading edges. 
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Figure 4. - indicia! aerodynamic characteristics of 
sinking wings with supersonic ieading edges. Mq-B. 
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0 .5 1.0 1.5 2.0 2.5 


Chord lengths traveled >T o 

(b) Pitching moment about 
leading edge or apex. 


Figure 4. - Continued. 
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(cj Location of center of pressure 
with reference to the teading 
edge or apex. 


Figure 4.— Concluded. 



Chord lengths traveled, 

(q) Lift. 

Figure 5 .- tndiclal aerodynamic characteristics of 

wings with supersonic edges pitching about lead- 
ing edge or apex. 2. 
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I 



(b) Pitching moment about 
leading edge or apex. 

Figure 5 - Continued. 
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Chord lengths traveled, 


( c) Location of center of pressure 
with reference to the leading 
edge or apex. 


i 


Figure 5 .- Concluded. 
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y^-(x^Mot)//9 
y ~ - m(x 


-y~(x^Mot)//3 
-y~m(x^ Mot) 

- x^ s M^ 


m- ctnA 




-y~mx*n/^m* 


-y=-mx*tv*m^ 


ly-’-mx-f I y-mx-tv*m^ 


Figure 7.- £ quations of tines pertinent to the anotysis 
of the triangular wing with subsonic leading edges. 
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Figure 8 - Regions used in the discussion of the 


aspect-ratio rectangular wing. 



(b) Pitching moment about apex. 
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Figure 9— Indicia! aerodynamic characteristics of sinking 
triangular wings with slender plan forms. mMo= i/8. 
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Chord lengths traveled, 


la) Lift. 



(b) Pitching moment about apex . 


Figure 10.^ Indicial aerodynamic characteristics of tri- 
angular wings with slender plan forms pitching about 

apex, mMo- i/8. 
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